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Abstract 

A hybrid numerical technique is presented for a characterization 
of the scattering and radiation properties of three-dimensional cavity 
arrays recessed in a ground plane. The technique combines the finite el- 
ement and boundary integral methods and invokes Floquet’s represen- 
tation to formulate a system of equations for the fields at the apertures 
and those inside the cavities. The system is solved via the conjugate 
gradient method in conjunction with the fast Fourier transform (FFT) 
thus achieving an O(N) storage requirement. By virtue of the finite el- 
ement method, the proposed technique is applicable to periodic arrays 
comprised of cavities having arbitrary shape and filled with inhomoge- 
neous dielectrics. Several numerical results are presented, along with 
new measured data, which demonstrate the validity, efficiency and ca- 
pability of the technique. 
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1 Introduction 


Recently, a hybrid finite element technique (FE-BI) was proposed for a 
characterization of the scattering and radiation properties of several three- 
dimensional cavity-backed structures including microstrip patch antennas 
and arrays [1]— [3] • The technique combines the finite element method with 
the boundary integral equation to formulate a system suitable for solution 
via the conjugate or biconjugate gradient method in conjunction with the 
fast Fourier transform (FFT). By virtue of the finite element method, the 
proposed technique is applicable to complex structures such as those involv- 
ing inhomogeneous dielectrics, conducting and resistive patches, feed probes 
and impedance loads. Accurate results have already been obtained for scat- 
tering and radiation by cavities, slots and microstrip patch antennas and 
these have demonstrated the method’s capability. 

In this paper, we develop the aforementioned finite element-boundary 
integral technique for scattering and radiation by infinite cavity arrays. Al- 
though, at least in principle, the technique is suitable for this application, 
there are new issues which must be addressed. For example, the applica- 
tion of Floquet’s theorem leads to a system which is non-symmetric. In 
addition, the matrix elements are dependent on the incidence angle in the 
case of scattering or the scan angle if the array is treated as a radiator. 
Consequently, for such a system an iterative solver is preferable over the 
usual direct solvers. Another issue addressed in this paper concerns the 
approximation of the finite array by a truncated infinite array. This kind 
of approximation has been adopted by necessity to make practical use of 
the infinite array analysis. However, its validity has not been examined for 
three-dimensional arrays and the two-dimensional studies are not conclu- 
sive [4]. Herewith we present three-dimensional comparisons of measured 
and calculated patterns which are perhaps the first to provide a measure of 
this approximation. 

The paper begins with a short description of the standard equivalence 
principle for subdividing the exterior and interior computational region. Us- 
ing Floquet’s theorem the exterior fields are formulated and expressed by 
an infinite sum involving the spectral representation of the exterior region’s 
Green’s function. This gives a set of discrete equations for the fields at the 
cavity interface. The interior fields are formulated via the finite element 
method leading to a sparse system of equations. By invoking tangential 
field continuity, this system is combined with the discrete set of equations 
developed for the exterior region. The final system is solved via the conju- 
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gate gradient method and certain restructuring is carried out to achieve an 
0(N) memory demand. Finally, a number of numerical computations are 
presented which demonstrate the capability and accuracy of the technique. 


2 Formulation 

The geometry under consideration is illustrated in Fig. 1, where a peri- 
odic array of cavities is recessed in an infinite ground plane. Each cavity 
is identical but may be of arbitrary shape and filled with inhomogeneous 
material. Also, the entire geometry may be covered with a thin dielectric 
layer. For scattering computations the excitation is assumed to be a plane 
wave and particular attention is given to the case of transverse electric (TE) 
and transverse magnetic (TM) incidences. For the radiation problem, the 
assumed excitation is a distributed current or current filament placed inside 
the cavities. 

A standard approach to formulate the scattering and radiation asso- 
ciated with the subject array is to employ the equivalence principle and 
close each cavity aperture with a perfect conductor. Above the aperture, 
an equivalent magnetic current M(x,y) = E(x,y) X z is introduced and by 
invoking tangential electric field continuity across the aperture, the appropri- 
ate equivalent sources for the fields in the cavity region axe — M(x,y) placed 
just below the aperture. To determine the equivalent currents M(x,y), the 
fields in each region are represented appropriately and by enforcing tangen- 
tial magnetic field continuity, a set of equations are derived for the solution 
of M(x, y). 

We begin the formulation by first invoking Floquet’s representation to 
express the magnetic current over the cavity apertures by the Fourier sum 

OO 

M(x, y) = (1) 

p,<7=-oo 

where 

M p , = JJ s M(x,y)i>* q ( x ,y)dxdy (2) 

■0p«j(x, y) = e , ' ( *" ,+fc *" v) . (3) 

In this, T x and T y denote the periodicities of the structure in the x and y 
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directions, respectively, S is the aperture of a unit cell and 

27 r 


— rp P “h k x , 

X 


, 27T - 

k yi = 7fr<l + k y 


with k' x = fc 0 sin 0* cos <f>' and k x y = k 0 sin d' sm <f>' . Also, {O' ,<f>') are the 
incidence or scan angles and k 0 = 2w/X is the free-space wavenumber. On 
invoking image theory, the field produced by M(x, y) in the upper half space 
can now be written as 


roc roc _ 

H sc (r) = -2jk 0 Y 0 / G{ r, r') • M(r') dS' 

J — OO J — oo 


( 4 ) 


where Gdenotes the pertinent dyadic Green’s function and Zo = 1/Fo is the 
free space intrinsic impedance. Assuming that G can be written in the form 

G = [°° r 8(*r, ky) dkx dky , ( 5 ) 

Alt 1 J — oo J — oo 

where f(z) is a function which reduces to 1 at z = 0, a substitution of this 
expression along with (1) into (4) yields 

oo ~ 

H sc (r) = -2jk 0 Y 0 £ G{k xp ,k yg ).M pq ^ x+k ^f{z). (6) 

p,(?=-oo 

Since H 5C (r) represents only the magnetic field due to the equivalent mag- 
netic sources on the ground plane, the total field in the upper half space is 
given by 

H(r) = H inc (r) + H ref (r) + H sc (r) (7) 

where H' nc is that generated by external sources (incident field) and H ref is 
the corresponding reflected field by the ground plane without the aperture. 
Both are nonzero only for scattering computations. 

We now consider the fields inside the cavity. In accordance with the 
variational principle, these must satisfy the variational equation [1]— [3] 

6F(E) = 0 (8) 

where the functional F is given by 

F(E) = \IIIv [^-( VxE )'( VxE )-^rE-E] dV 


+ f,h z " j 

-jk 0 Z 0 JJ 


int E - — M int -(V x E) j dV 

Mr 

M • H dS. 


( 9 ) 


4 



Here, V represents the cavity volume and (e r ,/i r ) denote the relative permit- 
tivity and permeability of the material filling the cavity. Also, (J ,nt ,M ,nt ) 
are the sources within the cavity and are nonzero only when the cavity 
houses an antenna configuration. The given functional is stationary about 
E and was derived on the assumption that the cavity’s aperture is closed 
with a perfect conductor and in the presence of the magnetic current sheet 
placed just below the closed aperture. 

To solve for the electric field via (8), it is necessary to replace the mag- 
netic field H in terms of the electric field or the magnetic current M = E x z. 
Since a knowledge of H is required only over the boundary/aperture sur- 
face 5, Maxwell’s equations cannot be used and this makes it necessary to 
introduce a different/independent condition relating E and H on S. Noting 
that tangential H must be continuous across the aperture, H in (9) can 
be replaced by the exterior magnetic field (7) for r € S. Doing so, the 
functional (9) can be written as 


F(E) = \JJJ v -J-(VxE)-(VxE) — *fc r E E dV 


2 

+ 


jk 0 Z 0 y nt • E M int • (V x E) 

VI Hr 

oo ~ 


dV 


— 2/$ // M • T G(k xp ,k yq )-M pq ^ x+k ^dS 

p,q=—oo 

— 2jkoZo Jj M • H inc dS. (10) 


Since M = E x z , this functional actually involves only the electric field and 
can therefore be discretized via the finite element method for a solution of 
the cavity and aperture fields via (8). 


3 Discretization and solution 


The finite element discretization of the first, second and last integrals of (10) 
are straightforward once the basis functions are chosen, as described in [2]. 
We will therefore concentrate on the discretization of the third integral 


rr °° - 

I = II M* ^ ] G{kxp,kyq) 

•'•'S p,q =- oo 


M pq e j V'*r x+k » 


dS. 


( 11 ) 


We first divide the region occupied by the periodic cell into small rectangles 
of dimension Ax X Ay. Within each subdivision, the equivalent current may 
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then be expanded using rooftop basis functions as 

M/2-1 N/2-1 

M x (x, y) = 12 M xmn A m+ i/ 2 {x)U n (y) (12) 

m——M/2 n=—N/2 
M/2-1 N/ 2-1 

My{x,y) — My mn II m (i) (13) 

m=—M/2 n=—N/2 


where 

f 1 — |x — mAx|/Ax, |x — mAx| < Ax 
Am(z) = | Qi | x _ mAx | > Al 

n I 1 ” mAx l < Al / 2 

m ' 1 \ 0, jx — raAx| > Ax/2. 


(14) 

(15) 


In the above, M and N denote the number of cells along the x and y 
directions, respectively, and M xmn and M ymn are the unknown constant 
coefficients of the expansion. Substituting these into (11), it is not difficult 
to show that / can be written as 


I — T x T y (iM xpq + yM ypg j -G (k xp , k yq ) • ^ iM X pq + yMypq'j 

P,Q 

where 

M/2-1 Nj 2-1 

M xmn e~i( 2 * p ' m / M+2 * q ' n / N } 

m= — M /2 n=— N/2 

- iWSM 


M, 


M, 


yvt 


MN 

M/2—1 N/ 2-1 

12 12 Mymnt 

m= — M /2 n= — N/2 


- j(2 j jrp , m/M+2nq , n/N ) 


(16) 


(17) 


(18) 


with p' = p + k x T x / 2 tt and q' = q + k'yT y /2TT. 

To construct the system implied by (8), it is necessary to differentiate F 
with respect to the field expansion coefficients and set the result to zero. For 
each coefficient, a single equation is obtained to make-up the overall system. 
The equations resulting upon differentiation with respect to the interior 
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cavity fields have been developed in [l]-[3]. As part of this process, when 
differentiating with respect to the aperture fields we encounter derivatives 
of I and from (16) these are found to be 


dl 

dM, 




dl 


dM, 


yst 


T x T\ 
(MN) 


rW G xx (k xp , k yq ) sine 4 (p) sine 2 fe) 

> p,i \ / \ 


x ^ ^f^ ne -jC2np'm/M+2*q'n/N) € H2xp's/M+2*q't/N) 


+ 


T X Ty 


(MN) 2 


^3 Gxy(kxp, k yq ) sine 3 (p) sine 3 (?p) q ' /N) 

P,Q 

X J3 M ymn e ~j( 2r P' m /M +2irq'n/N) gj(2irp's/M+2irq't/N) (19) 

(MN) 2 ^ Gyxik^ kyq) sine 3 ) sine 3 ) e Mp /M q /N) 

x ^ M ^ m ^ e -j(2Trp'm/M+2Tq'n/N) e j{2Tp's/M+2Trq't/N) 


+ 


T X Ty 


X3 G yy (k xp , kyq) sine 2 sine 4 (p) 

P.9 \ / \ / 

x ^ ^^^ e -j(27rp , m/Af+23T<7 , n/Ar) e j(2irp , 5/Af+27rg , i/A/’)^ (20) 


(A/jV ) 2 


m,n 


The resulting system can be solved on imposing the boundary condition at 
the cavity’s metallic walls which amounts to setting the tangential electric 
fields to zero at those walls. One of several solution algorithms can be used. 
However, since some of the matrix elements are in terms of the incidence or 
scan angles (see definitions of k XJ> and k yp ), direct solution methods such as 
the Gaussian elimination or LU decomposition methods are unattractive. 
This makes iterative solvers more suitable and we will specifically employ 
the conjugate gradient method. 

In implementing the conjugate gradient algorithm, one needs to compute 
matrix-vector products which involve the computation of (19) and (20). This 
can be efficiently accomplished via the FFT and without a need to generate 
the associated full submatrix. With this goal in mind, we rewrite (19) and 
(20) as 


dl 

dM xat 


ht 


TJTy 

(MN) 2 


M/2-1 N/2—1 

"y ' ^ ' A xx (k x p", kyq") 

p"=-M/2 q"=-N/2 
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M/2-1 N/2-1 

X £ £ 

m=— M/2 n=-N/2 

M/2-1 N/2—1 

+ fM TV ) 2 ^ ^ A x y( k xp ", k y t ]") 

> p"=— M/2 q"= — N/2 

M/2-l N/2-1 

X E E ^ymn/mne- J ' 2 ’ r(p " m/M+9 " n/N) (21) 

m=— M/2 n=— N/2 

„ „ M/2-1 N/2-1 

J-x-L-u V > V ^ T ✓* » \ 


dM~t fsi ~ (mV) 2 53 E 

aiK, y 4t ' ' p"=-M/2q"=-N/2 

M/2-1 N/2-1 

m=— M/2 n=— N/2 

M/2-1 N/2-1 

+ (Mm2 E E ^yy(^xp»’ fc y9") 

' > p"=— M/2 q"=—N/2 

M/2-l N/2-1 

x E E ^ymn/ m ne- j2 " (p " m/M+9 " n/N) (22) 

ra=— Af/2 n=—N/ 2 


where / mn = and 

oo 

? ^■yg // ) ~~ ^ ^ ^xx(^xp ;/ S sillC 




iv ; 


~ 00 - /p //; 7r\ 

A xy (k X p"i,k yq i") = E G X y(k X pin, kyqtii) sine 3 y ^ J s ‘ nc \J~n~ } 


X€ Mp"'/M-q"'/N) 


00 ~ / p //f 7T \ / ^ /// 7T \ 

A yx (k X piii , kyqiit} = E G yx (k X p,",kyqm) sine 3 sine 3 

U,tr=— OO ^ ^ ' 


xe -jx(p'"/A/- 9 "7N) 


°° /p /w 7 r\ fq"'x\ 

Ay y (k xp i" . kyqt'i) = E Gyy(k X pm, kyqi") sine J s * nc J 


with p"' = p" + uM and g"' = g" + t;jV. In accordance with the definition 
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of the discrete Fourier transformation, (21) and (22) can be symbolically 
written as 

= M N ^ 1 {A xr (k X p", kyq») ° T {M xmn f mTl } 

+ A xy (k xp ", k yq » ) o f{M ymri f mn }} (27) 

— ^ 4t } = \.Ay x (kxp"ikyq") ° T {M xmn f mn } 

+ Ayy(k X p" , kyg") O J- {Afymn/mn}^ (^8) 

where T and denote forward and inverse Fourier transformation and 
the symbol o implies the Hadamard product. This result demonstrates the 
convolutional form of the boundary integral and is essential in retaining an 
0( N) memory requirement for the solution algorithm. Because of (27) and 
(28), the generation of the full submatrix for the aperture fields is eliminated. 
Details relating to the implementation of the conjugate gradient algorithm 
in conjunction with the FFT operations required in (27) and (28) can be 
found in the literature [5]. 

4 Numerical results 

The formulation described above has been implemented and tested for sev- 
eral scattering and antenna geometries. In the following, we present some 
results which demonstrate the validity and capability of the method. By 
practical necessity, all of the computed scattering patterns will be for a fi- 
nite size array. They were obtained by multiplying the scattering pattern 
for a single unit cell of the infinite array with the standard array factor. 

First, we consider a 21 x 21 array of rectangular cavities of size 0.9A x 
0.4A x 0.1A. The cavities are 1A apart in the x direction and 0.5A apart in 
the y direction. The approximate backscatter RCS for this finite array is 
given in Fig. 2 and is compared with results obtained via the so-called mode 
matching technique [6]. Both RCS patterns are virtually the same, but 
this only validates the infinite array solution since the finite array pattern 
computed via the mode matching technique was also obtained from a similar 
approximation of the infinite array solution. 

To validate the approximate finite array solution, it is necessary to com- 
pare our results with an exact solution which is based directly on the finite 
array model. No such results are available and thus for this purpose, a 5 x 5 
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array was manufactured consisting of 3 cm x 3 cm x 0.45 cm cavities placed 
4 cm apart from each other in both the x and y directions. This array was 
mounted on a diamond-shaped low cross section body and the backscatter 
RCS was measured as a function of incidence angle in the principal plane 
perpendicular to the structure for both Eg and Ej, polarizations. Measure- 
ments were collected at several Af-band frequencies and in all cases a good 
agreement was observed between theory and experiment. A typical compar- 
ison is given in Fig. 3 at 9.1 GHz and we remark that no measured data are 
displayed near broadside and grazing due to the obvious interference from 
the test platform. 

In addition to scattering computations, the proposed method can also 
be employed for the analysis of a periodic microstrip patch array with each 
radiating element situated in the cavity. The microstrip patch array in the 
absence of cavities, i.e., on an infinite substrate, has, of course, been studied 
extensively [7], [8]. However, it has been proposed [9] to place metal walls 
between the radiating elements to reduce mutual coupling which in turn 
could cause severe mismatch and blindness. Aberle [10] applied this idea 
to microstrip patch arrays using a modal technique, specifically developed 
for a metal circular cavity and uniform filling. By virtue of the finite ele- 
ment method, the presented technique allows much more flexibility in cavity 
shape and filling. While an extensive investigation of this problem is beyond 
the scope of this paper, here we present two examples which illustrate the 
method’s capability. For this purpose, let us consider the array of cavity- 
backed patches illustrated in Fig. 4. In the first example, we assume that 
a 0.3A x 0.3A rectangular patch is placed at the aperture center of a cavity 
0.45A x 0.45A x 0.02A in size and filled with a dielectric having e r = 2.8 
(the periodicity of the cavity array is one half of a wavelength in either 
directions). The computed active input resistance and scanning reflection 
coefficient for this array are displayed in Fig. 5 as a function of scan angle. In 
the second example, we assume that a 1.8 cm x 1.8 cm rectangular patch is 
placed at the aperture center of a cavity 3.3 cm x 3.3 cm x 0.318 cm in size 
and filled with a dielectric having e r = 2.33 (the array’s periodicity is now 
3.6 cm in either direction). The calculated scanning reflection coefficient for 
this array at 4.9 GHz is shown in Fig. 6. It is clear from the pattern in Fig. 6 
that nearly total scan blindness occurs at 6 = 43.8° in the E-plane but no 
blindness is observed in the H-plane. Note that the blindness angle can be 
predetermined from grating theory and for the uncoated array it occurs only 
when T x or T y is greater than one half of a wavelength. When the array is 
coated with a dielectric layer, additional blindness angles may occur which 
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can be also predetermined from surface wave theory [7]. 

5 Conclusions 

In this paper, we presented a hybrid numerical method for scattering and 
radiation analysis of a three-dimensional cavity array recessed in a ground 
plane. The solution method combined the finite element and the boundary 
integral formulations along with Floquet’s field representation to generate 
a matrix system for the cavity and aperture fields in a periodic cell. This 
system was solved via the conjugate gradient method in conjuction with 
the FFT for an efficient evaluation of the required matrix-vector products. 
Numerical results were presented to demonstrate the validity, accuracy and 
capability of the method and these included backscatter RCS patterns, input 
impedance and scanning reflection coefficient computations for microstrip 
patch arrays. 


Appendix: Dyadic Green’s functions 


To implement the formulation described in the paper, we need explicit ex- 
pressions for the dyadic Green’s function in the spectral domain. When the 


ground plane is uncoated, these are given by 

G xx (k x ,k y ) = (&o ~ kl)/(2jk z ) (29) 

Gyy(k x ,ky) = (l£ ~ fy / (2 jk z ) (30) 

G xy (k x ,ky) = G xy (k x , ky) = -k x ky/(2jk z ) (31) 

where k z = - k^. For a coated ground plane, the required expres- 

sions are 

G xx (k x ,k y ) = {erklk^ - k 2 zl klT 2 ) /(2k 2 0 k zl p 2 ) (32) 

G yy (k x ,ky) = (e r k 2 klT l -k 2 zl k 2 y T 2 )/(2k 2 k zl p 2 ) (33) 


G xy {k x , ky) = G yx = -k x ky (/ k 2 zl T 2 + e r k 2 T 1 ) /(2 k 2 k zl (3 2 ) (34) 


where 


€ r kz 2 sm k z \d — jk z i cos k z \d 
e T k z 2 cos k z \d + jk z \ sin k z \d 
k z2 cos k z id + jk z \ sin k z \d 
k Z 2 sin k z \d — jk z i cos k z \d 


(35) 

(36) 
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k ii = e r ko - /? 2 , = ^0 - /? 2 > = k l + and d is the thickness of the 

dielectric layer. 
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Figure captions 

Fig. 1. Problem geometry. 

Fig. 2. Backscatter RCS oqq of a 21x21 array of rectangular cavities in a 
ground plane. For this computation T x = lA,T y = 0.5A, the cavity 
size is 0.9A x 0.4A x 0.1A, and <j> — 20°. 

Fig. 3. Principal backscatter RCS pattern of the 5 x 5 cavity array. For this 
computation T x = T y - 4 cm, the cavity size is 3 cm x 3 cm x 0.45 
cm and /=9.lGHz. (a) E$ pattern, (b) pattern. 

Fig. 4. Array of patches residing on the surface of individual rectangular cav- 
ities. 

Fig. 5. Scanning reflection coefficient and active resistance as a function of 
scan angle for an infinite cavity-backed patch array. For this compu- 
tation T x = T y = 0.5A, the cavity size is 0.45A x 0.45A x 0.02A and the 
patch size is 0.3A x 0.3A. The cavity is filled with a dielectric having 
£ r = 2.8. 

Fig. 6. Scanning reflection coefficient as a function of scan angle for an infinite 
cavity-backed patch array. For this computation T x = T y = 3.6 cm, 
the cavity size is 3.3 cm x 3.3 cm x 0.318 cm and the patch size is 1.8 
cm x 1.8 cm. The cavity is filled with a dielectric having e r = 2.33. 
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a/r (dB) 


Backscatter RCS of a 21x21 Cavity Array 



0 (degrees) 


Fig. 2 Backscatter RCS of a 21x21 array of rectangular cavities in a ground plane. 
For this computation T x = 1 A,T y = 0.5A, the cavity size is 0.9A x 0.4A x O.lA, and 
<f> = 20°. 
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Fig. 3 Principal backscatter RCS pattern of the 5x5 cavity array. For this computation 
T x = T y = 4 cm, the cavity size is 3 cm x 3 cm x 0.45 cm and /=9.lGHz. (a) Eg 
pattern, (b) E 4 , pattern. 18 






Top view 
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Cross-sectional view 


Fig. 4 Array of patches residing on the surface of individual rectangular cavities. 
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(b) Scanning angle (degrees) 

Fig. 5 Scanning reflection coefficient and active resistance as a function of scan angle 
for an infinite cavity-backed patch array. For this computation T x = T y = 0.5A, 
the cavity size is 0.45A x 0.45A x 0.02A and the patch size is 0.3A x 0.3A. The cavity 
is filled with a dielectric having e r == 2.8. 
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Fig. 6 Scanning reflection coefficient as a function of scan angle for an infinite cavity- 
backed patch array. For this computation T x = T y = 3.6 cm, the cavity size is 3.3 
cm x 3.3 cm x 0.318 cm and the patch size is 1.8 cm x 1.8 cm. The cavity is filled 
with a dielectric having e T = 2.33. 
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